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HIGH ALTITUDE VENUS OPERATIONS CONCEPT TRAJECTORY 
DESIGN, MODELING, AND SIMULATION 
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A trajectory design and analysis that describes aerocapture, entry, descent, and 
inflation of manned and unmanned High Altitude Venus Operation Concept 
(HAVOC) lighter-than-air missions is presented. Mission motivation, concept of 
operations, and notional entry vehicle designs are presented. The initial trajecto- 
ry design space is analyzed and discussed before investigating specific trajecto- 
ries that are deemed representative of a feasible Venus mission. Under the pro- 
ject assumptions, while the high-mass crewed mission will require further re- 
search into aerodynamic decelerator technology, it was determined that the un- 
manned robotic mission is feasible using current technology. 

INTRODUCTION 

For decades, Mars has been the prime destination for a variety of robotic science missions as 
well as a target for human exploration. However, there are a number of challenges in executing a 
manned Mars mission: transit times are significant (as much as a year each way), the surface en- 
vironment is harsh with low ambient temperatures and pressures, and the surface radiation envi- 
ronment is hazardous to crew health. Venus, conversely, has a shorter transit time (approximately 
100 days outbound and up to 300 days inbound) with comparable requirements on in-space pro- 
pulsion systems. While the surface of Venus presents an environment with temperatures and pres- 
sures too great to permit surface platforms to survive for any length of time with current technol- 
ogies, the atmosphere at approximately 50 km altitude is one of the more hospitable environments 
in the solar system; the ambient temperature and pressure are much closer to that on the surface 
of Earth. Radiation protection for the crew is provided by the atmosphere, and a life-supporting 
mixture of gas similar to air could dually serve as a lifting gas. The Venus Exploration Analysis 
Group^ (VEXAG) and the Planetary Science Decadal Survey^ have also identified planetary sci- 
ence questions that may only be answered from operations in the Venusian atmosphere. 

The present work is in support of a NASA Langley Research Center (LaRC) study to consider 
a conceptual design of an atmospheric platform that would be capable of supporting both human 
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and robotic missions. This study of a High Altitude Venus Operational Concept (HAVOC) con- 
siders the goals, requirements, and systems to explore Venus using both robotic precursor and 
manned vehicles. A lighter-than-air (LTA) vehicle, or airship, offers the capability of meeting the 
study's objectives. Analysis of the aerocapture and entry, descent, and inflation trajectories is crit- 
ical to developing the concept for the airship and its associated atmospheric entry systems. 

The present paper will focus on the design, modeling, simulation, and analysis of the two pri- 
mary segments of the HAVOC trajectories: aerocapture, and entry, descent, and inflation (EDI). 
During aerocapture the vehicle takes advantage of the drag created by flying through the atmos- 
phere to reduce the speed of the vehicle, without the use of a large propulsive maneuver. Bank 
angle modulation is used to adjust the lift vector and guide the vehicle toward the desired target 
apoapsis, using a small propulsive burn to clean up any targeting errors. Once the vehicle is in 
the desired circular orbit, atmospheric entry may begin. During entry, the entry vehicle (EV) is 
guided through the middle and lower atmosphere and may be maneuvered to prevent excessive 
heating and aerodynamic loads. Once the EV has slowed to supersonic velocities, the descent 
phase begins when an aerodynamic decelerator, such as a parachute or ballute, is deployed to fur- 
ther reduce descent rate. The inflation phase begins when the EV velocity is sufficiently reduced 
such that the airship may be exposed to the oncoming flow. The airship is pressurized during par- 
achute descent and reaches a nominal altitude of approximately 50 km. Figure 1 shows the 
HAVOC concept of operations from aerocapture to inflation and operations. Though not shown 
in the figure, the manned mission requires two vehicles to undergo aerocapture: the crew and air- 
ship/cargo are launched separately and arrive at Venus in separate vehicles, then rendezvous in 
orbit before performing the EDI maneuver. 
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Figure 1. HAVOC Concept of Operations. 

TRAJECTORY DESIGN, MODELING, AND SIMULATION 

Design and analyses of the entry and descent phases of the trajectories were conducted using 
the Program to Optimize Simulated Trajectories II (POST2).^ '^ Three-degree-of-freedom simula- 
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tions run with POST2 used aerodynamic and aeroheating databases generated with CBAERO for 
the scaled mid-L/D “ellipsled” from the NASA Entry, Descent, and Landing Systems Analysis 
(EDLSA) study.^ 

The nominal atmosphere was modeled using VenusGRAM 2005^ assuming no winds. The 
Venus gravitational field was modeled up to the J4 zonal harmonic. Entry interface was defined 
as a geodetic altitude of 200 km, and the entry point was defined as the intersection of the equator 
and the Venus prime meridian (0° latitude, 0° longitude). 

Figure 2 shows nominal atmospheric density, pressure, and temperature profiles of Venus ob- 
tained from VenusGRAM 2005. The conditions at 50 km altitude are roughly equivalent (though 
warmer) to atmospheric conditions at sea level on Earth; the nominal density, pressure, and tem- 
perature are 1.5948 kg/m^, 106,679 Pa (1.0528 atm), and 350 K (170°F), respectively. Airship 
operations are therefore designed for 50 km. 



Figure 2. Nominal atmospheric profile of Venus, mid to high altitudes. 


The sharp density gradient, coupled with the very high entry mass of the manned mission, 
present significant trajectory design challenges. Among these are: 

1. Atmosphere skip-out. Skip-out results when the EV is traveling at too high a velocity (or 
has too much lift) to be “captured” and exits the atmosphere at a reduced speed, and may 
or may not re-enter. When controlled, this may be used as a technique for aerocapture 
and aerobraking. 

2. Excessive lofting. Lofting is a repeated increase and decrease in altitude due to excessive 
lift that may result in atmosphere skip-out, and is usually more pronounced in trajectories 
with shallow entry flight path angles. 

3. Excessive aerodynamic forces (g-loads). Excessive g-loads may injure crewmembers and 
can be fatal at high levels over extended periods. For this reason, NASA has established 
g-load limits for human crews that are functions of length of exposure and orientation of 
the crew relative to the acceleration vector.^ 

4. Excessive aerodynamic heating. Atmospheric entries are characterized by high heat rates 
and loads that require the use of thermal protection systems (TPS) to dissipate heat and 
protect the EV. It is advantageous to minimize heat rates and loads, thereby reducing 
structural mass. 
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To address these ehallenges, a variety of different trajeetories were analyzed by establishing a 
parameter design spaee and varying speeifie independent parameters relating to the trajeetory. 
The result was thousands of simulated trajeetories from whieh parameters of interest were ex- 
traeted and analyzed for trends. The independent parameters of interest were entry flight path an- 
gle, total angle of attaek, and bank angle. Dependent parameters of interest ineluded peak g-load, 
peak heat rate, and total heat load. A trajeetory was eonsidered “valid” if it did not skip out and if 
the peak g-load did not exeeed predetermined limits. Peak heat rates were monitored, though not 
used as a seleetion eriteria, sinee the thermal proteetion system is designed around these values. 

Initial analyses were eondueted with an unguided vehiele with a theoretieal eontrol system that 
maintained a eonstant total angle of attaek and bank angle through the trajeetory. A target “land- 
ing” loeation was not eonsidered due to the highly mobile nature of the payload/airship, whieh 
does not require injeetion at a preeise altitude or latitude. Later analyses implemented a simple 
logie algorithm to modulate bank angle to mitigate aerodynamie loads. 

Aeroshell and Airship Modeling 

Initial trajeetory designs assumed diserete aeroshells for both aeroeapture and EDI, that is, the 
aeroeapture aeroshell would be jettisoned prior to EDI. As the design of the payload matured and 
the overall mass of the system inereased, it was determined that the initial design would not be 
eapable of supporting the heating values predieted by the simulation for either aeroeapture or 
EDI. Other deeelerator eoneepts were explored and a rigid mid-L/D vehiele was seleeted that 
would be able to aeeommodate the mass of the system as well as the two heat pulses that result 
from two distinet passes through the Venus atmosphere. As will be diseussed in the present work, 
the eurrent designs feature a single aeroshell for both aeroeapture and entry with a thermal protee- 
tion system tailored to the modeled surfaee heat distribution. 

Figure 3 shows notional aeroshell and airship designs. Relevant design parameters for the en- 
try vehiele aeroshells are listed in Table 1. 



Figure 3. Notional entry vehicle aeroshell (left) and airship (right) designs. Images are not 

to relative scale. 


Table 1. Entry vehicle parameters. 


Parameter 

Unit 

Unmanned 

Manned 

Diameter 

m 

4.7 

10.0 

Entry mass 

mt 

4.537 

133.9 


Both HAVOC airship concepts (unmanned and manned) were sized using the SolFlyte-LTA 
sizing utility. SolFlyte was developed within the Aeronautics Systems Analysis Branch (ASAB) 
at NASA Langley Research Center for modeling vehicles and extended-duration missions carried 
out utilizing solar-electric LTA concepts. The model provides robust, general, rapid sizing esti- 
mates and iterative parametric design capability. Key input and output specifications from the 
LTA sizing utility were used as the basis to formulate the HAVOC inflation and descent model- 
ing capability. 
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The existing Earth-based airship sizing and operational methodology was slightly modified to 
refleet the speeifie requirements of the Venus-based HAVOC eoneept. Determining airship 
weight and volume required implementing eonversion faetors related to differenees in gravity, 
atmospherie moleeular weight, and buoyaney. The speeifie buoyaney of the lifting gas was used 
to transform the Earth-based design atmosphere density of 1.15 kg/m^ to a Venus -based density 
of 1.0914 kg/m^ that eorresponds to altitudes near 50 km. Caleulated design volumes of the 
HAVOC airship eoneepts were 1,118 m^ and 77,521 m^ for the robotie and human missions, re- 
speetively. Caleulated total masses of the airships, ineluding lifting gas and tanks were 1,382 kg 
and 95,776 kg for robotie and human missions, respeetively. 

A signifieant design eonsideration beeame apparent during initial assessments of the human 
mission. Airship designs do not typieally inelude a requirement to transport lifting gas from one 
planet to another. In faet, ground-based inflation faetors into airship design as an operational eon- 
straint. For the HAVOC mission, however, transporting the lifting gas mass to Venus presents a 
signifieant payload eonstraint. The HAVOC erewed mission requires transporting 8,200 kg of 
helium lifting gas, not ineluding tanks. If it were desirable to use a mixture of gases similar to air 
as both the lifting gas and a means of life support, 59,400 kg of the gas mixture would be re- 
quired. Thus, for a human mission, redueing the spaeeeraft payload mass dietates utilizing the 
lowest moleeular weight lifting gas eapable of life support. Though not eonsidered in the present 
analysis, it will be of interest to investigate the use of a dual-purpose helium-oxygen mixture for 
the manned mission. 

Aerodynamics and Aerothermodynamics 

The vehiele outer mold line (OML) eonsidered for the aeroeapture and entry phases of flight is 
an axisymmetrie ellipsled based on that of the ELDSA^ and Design Referenee Arehiteeture 5.0 
(DRA-5)^ studies. The eonfiguration may be deseribed as a right eireular eylinder with a hemi- 
spherieal nose eap and flat base, where the baseline exploration-elass or manned vehiele has a 
diameter of 10 meters and overall length of 30 meters, and the robotie-elass vehiele dimensions 
are downsealed to 47% of the baseline. Note that the aeroshell design in Figure 3 is notional and 
shown for illustrative purposes, and is more eomplex than the aerodynamie model used in the 
analysis deseribed here. 



Figure 4. Mach number and dynamic pressure combinations used for CBAERO solutions. 
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Aerodynamic and aerothermodynamic databases for both exploration- and robotic-class ellip- 
sleds were generated using the Configuration Based Aerodynamics (CBAERO) code.^ Solutions 
were computed over a matrix of Mach number and dynamic pressure combinations ranging from 
free-molecular aerocapture to low supersonic continuum regime conditions as shown in Figure 4, 
where the transition between free molecular and continuum computations was blended automati- 
cally based upon an internal calculation of Knudsen number. At each point, solutions were com- 
puted for total angles of attack ranging from 0° to 60° in increments of 10°. A Venus atmosphere 
model extending to 250 km altitude was used for atmospheric properties including temperature, 
pressure, ratio of specific heats, and molecular weight. Radiative heating environments were pre- 
dicted with a Tauber-Sutton correlation,^^ while fully turbulent flow and a fully catalytic wall 
boundary condition were assumed for conservative convective heating estimates. No control sur- 
faces were modeled, as the corresponding sizing and requisite control algorithms were considered 
beyond the scope of this study. Point checks for aeroheating and aerodynamics were computed at 
representative Mach 24 entry conditions using the LAURA^^ and OVERFLOW^^ CFD codes, 
respectively, and were in good agreement with CBAERO results. 

Aerocapture 

For the purposes of this study, the aerocapture maneuver is designed to take each vehicle from 
arrival speeds of approximately 10 to 12 km/s to the roughly 7 km/s velocity necessary to hit the 
desired orbit apogee. During the pass through the atmosphere, bank angle modulation is used to 
adjust the lift vector and ensure the vehicle hits the target orbit with minimal error. A simplified 
bank angle profile is implemented to emulate the flight performance of the vehicle while allowing 
for the rapid generation of trajectories. 

In order to understand the conditions experienced by each vehicle, a design space was devel- 
oped exploring combinations of input parameters including arrival velocity, target orbit, system 
mass, and lift-to-drag ratio (L/D). One design space was created for each of the three systems in- 
cluded in this study: one for the robotic mission and two for the manned mission (recall that the 
crew and airship/cargo arrive at Venus in separate vehicles). To fully characterize the design 
space, trajectories were generated which would produce minimum and maximum values of all the 
desired output parameters for each combination of inputs. This was completed by examining the 
aerocapture corridor. For each combination of inputs, the corridor is bound by flying one trajecto- 
ry with the lift vector pointing up the entire pass (bank angle of 0°) and then iterating on the entry 
flight path angle value which allows that trajectory to hit the apoapsis target, then repeating the 
process for a pass with the lift vector pointing down (bank angle of 180°). These are referred to as 
the steep and shallow trajectories, respectively, and bound the possible range of flight conditions 
experienced by the vehicle during aerocapture. One corridor was generated for each combination 
of inputs described above and this process was repeated for each of the three systems considered 
for this study. 

Trajectory outputs from the three sets of corridors were used to determine which set of inputs, 
for each vehicle, would produce the desired nominal conditions. Contour plots, demonstrated in 
Figure 5, were generated for each system to examine conditions such as peak g-loads, peak heat 
rates, and dynamic pressures. Using the mission constraints imposed on these conditions, one 
nominal trajectory was generated for each of the three systems to be used for TPS sizing. A sim- 
plified bank angle profile was implemented for the nominal trajectory, which involves flying the 
vehicle with the lift vector pointing directly up (bank angle of 0°) in order to bleed energy from 
the trajectory without incurring excessive g-loads and heating rates, and then reorienting the vehi- 
cle with lift vector pointing down (bank angle of 1 80°) once peak dynamic pressure has passed, in 
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order to continue to bleed energy until the desired exit velocity is reached. This method produces 
three nominal trajectories, one for each system, from which iterations on the design could begin. 
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Figure 5. Sample contour of output parameters for HAVOC robotic mission. 


Entry and Descent 

The entry and descent (ED) portion of the trajectory begins at entry interface (200 km) and, 
for the purposes of the present analysis, terminates at supersonic parachute deploy (Mach 2.1). 
The following assumptions were made to facilitate fast analysis of the ED trajectory design space: 

1 . A constant total angle of attack, selected to produce a desired amount of lift, is main- 
tained from entry interface to parachute deploy. Note that the EDLSA ellipsled produces 
maximum lift at a 55° total angle of attack.^ Angles of attack greater than this were ex- 
plored in the design space but were considered invalid. 

2. The entry vehicle has sufficient control authority to maintain such an orientation and 
flight path angle. 

3. A skip-out margin of 3° was specified to ensure skip-out would not occur. That is, a valid 
trajectory must have an entry flight path angle at least 3° greater than the flight path angle 
that results in skip-out. 

4. Convective and radiative heat margins of 1.30 and 1.50, respectively, were applied. 

5. Bank angle modulations are modeled as instantaneous step functions. 

An initial design parameter space was established to explore a range of trajectories bounded 
by lift-up and lift-down orientations. Velocities at entry interface were fixed at 7.2 km/s. Selected 
entry flight path angles ranged between -7.0° and 0.0°, total angles of attack ranged between 0° 
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and 60°, and bank angles ranged between 0° and 180°, where 0° is lift -up and 180° is lift-down. 
Eaeh trajeetory was run in POST2 and analyzed in MATLAB. Though dependent on the magni- 
tude of the lift veetor, whieh is in part determined by the total angle of attaek, lift-up trajeetories 
are generally eharaeterized by low g-loads, low heat rates, high total heat loads, and greater like- 
lihood of atmosphere skip-out. Lift-down trajeetories, eonversely, are generally eharaeterized by 
high g-loads, high heat rates, low total heat loads, and lower likelihood of atmosphere skip-out. 


Figure 6 and Figure 7 show peak g-load and peak heat rate variations aeross bank angles and 
total angles of attaek at a given entry flight path angle for the robotie and manned missions. It ean 
be seen that as the entry flight path angle beeomes steeper, the peak g-load at a given bank angle 
and total angle of attaek inereases. In faet, at an entry flight path angle of -7.0°, there is no eom- 
bination of bank angle and total angle of attaek that produees a peak g-load below 10 for either 
robotie or manned missions. Similarly, steeper flight path angles produee higher heat rates. In 
general, the “sweet spots” eorrespond to regions of shallow flight path angles with mid to high 
angles of attaek and near-lift-up orientations. 


Max g load variation (g) 
Entry flight path angle > -2 deg 



Total angle of attack (deg) 


Max g load variation (g) 
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Max g load variation (g) 
Entry flight path angle > -7 deg 



Total angle of attack (deg) 
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Figure 6. Robotic mission trajectory design space at selected flight path angles. Note that 
any values that exceed the range of the legend remain dark red or dark blue. 
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Figure 7. Manned mission trajectory design space at selected flight path angles. Note that 
any values that exceed the range of the legend remain dark red or dark blue. 


Many of the valid trajectories obtained from the initial design space, exhibited two character- 
istics that necessitated adjustment of the entry profile: lofting and excessive g-loads. To counter 
the lofting behavior, a basic control logic algorithm was implemented in which the vehicle bank 
angle was modulated to near 90° if the altitude began to increase, which effectively rotates the lift 
vector away from near-vertical, or lift-up, thereby reducing the increase in altitude due to lift. Ex- 
cessive g-loads were also mitigated by bank angle modulation, though in the opposite sense. That 
is, the bank angle was modulated to near 0° if the g-loads became excessive, which effectively 
rotates the lift vector towards the lift-up orientation. The result is a bank angle modulation profile 
that mitigates both lofting and excessive aerodynamic forces. To enable a clean deployment of 
the parachute, the final modulation is to a bank angle of 0° just prior to parachute deploy. Total 
angle of attack is also set to 0° to eliminate any lift. 

A disc-gap-band (DGB) supersonic parachute, which was the assumed aerodynamic decelera- 
tor of choice for the HAVOC missions, may in general be deployed at Mach numbers between 
2.25 and 1.00 and dynamic pressures between 900 and 200 Pa. The region in altitude -velocity 
space bounded by these limits is often referred to as the “Mach-q box” and is shown in Figure 8. 
Exceeding these conditions increases the probability that the parachute will either be destroyed 
due to excessive loading or exhibit areal oscillations that decrease efficiency of the parachute. 
Note the relatively low range of altitudes and velocities in the Mach-q box, the significance of 
which will be discussed later. 
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Figure 8. Venus Mach-q parachute box for (VenusGRAM 2005 atmosphere). 

The purpose of the entry and deseent portion of the trajeetory is to plaee the EV at a eondition 
that permits deployment of the supersonie paraehute, after whieh the airship may begin inflation 
under terminal paraehute deseent. To this end, the valid trajeetories were evaluated at Maeh 2.1 to 
establish if paraehute deployment eonditions were met. 

Airship Inflation and Terminal Descent 

The terminal deseent portion of the trajeetory begins after paraehute deploy. A HAVOC Ter- 
minal Deseent Model (HAVOC-TDM) was developed for analyzing how aerodynamie, buoyan- 
ey, and inertial forees eombine to fix the terminal veloeity of the vehiele eoneept during unpow- 
ered deseent. Existing trajeetory models were utilized for modeling the initial deseent from orbit 
to paraehute deployment. However, to ensure a smooth transition between separate modeling 
tools, HAVOC-TDM does provide the option of moderate-fidelity modeling of the aerodynamie 
eonditions prior to the time of paraehute deployment. Table 2 lists relevant assumed paraehute 
and airship parameters. 


Table 2. Parachute and airship parameters. 


Parameter 

Unit 

Unmanned 

Manned 

Supersonic Parachute 

Diameter 

m 

10 

24 

Supersonic drag coeff 

- 

0.9 

0.9 

Transonic drag coeff 

— 

0.85 

0.85 

Subsonic drag coeff 

— 

1.15 

1.15 

Airship 

Mass 

kg 

1,382 

95,776 

Packed volume 

m3 

20 

1258 

Inflated volume 

m3 

250 

77521 

Inflation rate 

m3/s 

1 

400 
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HAVOC-TDM determines the aerodynamie and buoyaney forees aeting on the transforming 
vehiele eonfiguration as a funetion of time at altitude inerements of 500 m. Table 3 is an example 
output from the model, ineluding the altitude, time segment, and eumulative deseent time. In ad- 
dition, a detailed table of the faetors required to ealeulate drag foree eomponents, and the buoy- 
aney, weight, veloeity, and aeeeleration is provided upon eompletion of eaeh analysis. 


Table 3. Sample inflation and terminal descent parameters. 



entry 

deploy 

descent 

descent 

descent 

descent 

End 1 

End 2 

End 


Entry 

Start 

Parachute 

Deploy 

Parachute 

Open 

Shell 

Cut 

Inflate 

Start 

Parachute 

Cut 

Inflation 
A Time 

Velocity 

zero 

DESCENT 

altitude, km 

69.5 

69.0 

68.0 

64.0 

64.0 

54.5 

52.0 

52.0 

52.0 

dtime, s 

58.55 

2.05 

5.48 

30.46 

35.94 

151.89 

190.7 

257.0 

257.0 

timeline, s 

58.55 

60.60 

66.08 

96.54 

96.54 

248.43 

287.2 

353.5 

353.5 


Various mission, vehiele, and subsystem input parameters are used to dietate the sequenee of 
operational events and to speeify the eonditions assoeiated with the onset of eaeh event. During 
eaeh step, the deseent eonditions are assessed to determine whieh events have already oeeurred 
and if another operational event should be initiated. Distinet modeling and ealeulations are ap- 
plied to eaeh phase of the deseent as a result of eharaeteristieally different vehiele eonfigurations, 
weights, and buoyaney forees. Thus, the deseent analysis model follows a pattern of 1.) identify- 
ing the deseent phase, 2.) ealeulating eonditions, 3.) assessing event onset, and 4.) ealeulating the 
time to deseend 500m. 

Onee the vehiele inflation phase begins, it is assumed that all entry shields and eapsules are 
jettisoned to reduee the total system weight, and that helium gas inflates the HAVOC eoneept at a 
speeified and eonstant volume rate. The buoyaney foree is ealeulated assuming an atmosphere 
moleeular weight of 43.58 g/mole (approximately 97% CO 2 and 3% N 2 ) and a helium lifting gas 
weight of 4.0 g/mole. Due to the signifieant mass of the lifting gas tanks, the lifting gas was as- 
sumed to be equally eontained in multiple tanks whieh would operate in sueeession and be jetti- 
soned upon depletion. This design assumption was determined to signifieantly slow the rate of 
deseent. 

Changing drag forees on the vehiele as a result of inflation are ealeulated, ineluding parasite 
and air-mass drag. The deseent drag, paraehute drag, buoyaney, and vehiele weight are used to 
determine the time required to deseend 500m for the existing eonditions. Updates to these faetors 
eontinue until the buoyaney to paraehute drag ratio exeeeds 90%, at whieh time the paraehute is 
released. The vehiele deseent rate inereases for a short time period after the paraehute is released, 
but the atmosphere density inereases to provide buoyaney inereases that exeeed those eaused by 
the additional lifting gas volume released during that same period. 

The deseent time to the point of neutral buoyaney is determined at the point that the buoyaney 
foree equals the vehiele weight. However, deseent may eontinue as a result of the momentum of 
the vehiele, whieh may or may not be inflated at that altitude. In this ease, as the vehiele reverses 
to aseend toward the higher altitude of equilibrated neutral buoyaney, the minimum altitude and 
the deseent time to that altitude are output. 

It is also possible that the inflation rate permits inflation at altitudes below the level of neutral 
buoyaney. In this ease, as the vehiele reverses to aseend toward the higher altitude of equilibrated 
neutral buoyaney, the minimum altitude and the deseent time to that altitude are output. It is also 
possible that the inflation rate permits inflation at altitudes somewhat above the level of neutral 
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buoyancy, in which case the inflated airship is assumed to safely descend to the neutral buoyancy 
altitude. 

Thermal Protection System 

The thermal protection system (TPS) was sized based on the nominal trajectory heating pro- 
files with a 25% margin added to the final mass. For the robotic mission, given a 4.7 m diameter 
cylinder and half-sphere nose cap with a total length of 14.1 m, and assuming that the TPS covers 
all of the spherical nose cap and half of the cylindrical body, the total area covered by the TPS is 
121.45 no?. For the manned mission, given a 10.0 m diameter cylinder and half-sphere nose cap 
with a total length of 30.0 m, and assuming that the TPS covers all of the spherical nose cap and 
half of the cylindrical body, the total area covered by the TPS is 549.78 m^. 

Two TPS candidates were investigated for both missions. The first, Heatshield for Extreme 
Entry Environment Technology (HEEET), is a dual-layer material, with a high-density outer “re- 
cession” layer woven in to a lower density “insulation” layer and 3D woven carbon fibers infused 
with phenolic resin. The HEEET heatshield thickness was calculated to keep the TPS -structural 
interface temperature below 315°C. HEEET is currently at a low technology readiness level 
(TRL), but technology development effort is underway to prepare HEEET for mission infusion. It 
is expected that HEEET will be robust enough to handle a dual heat pulse with a cold soak be- 
tween the pulses, which would be in the HAVOC mission profiles (aerocapture - cold soak - at- 
mospheric entry). 

The second TPS candidate was Phenolic Impregnated Carbon Ablator (PICA), which is a 
monolithic resin-infiised fiberform insulation. PICA requires a strain isolation pad (SIP) to iso- 
late it from the support structure. A thickness margin of 10% was added for recession uncertainty. 
The PICA heatshield thickness was calculated to keep the TPS -structural interface temperature 
below 250°C. Note that it is currently unknown if PICA can withstand a dual heat pulse with a 
cold soak between the pulses. 

For the TPS support structure, a 5.08 cm, 4.4 pcf aluminum honeycomb truss sandwiched be- 
tween 0.1 cm graphite/BMI facesheets was assumed for the robotic mission. The manned mission 
used a more robust 7.62 cm, 6.0 pcf titanium honeycomb truss sandwiched between 0.2 cm 
graphite/BMI facesheets. 


Section 4 



Section 1 Section 2 Section 3 

Q = Qmin @ X = 0 Q = Q = 0,^3, 


Figure 9. Moderate TPS tailoring. 

Three different sizing techniques were used: constant TPS thickness, moderate thickness tai- 
loring, and aggressive thickness tailoring. For the first, the TPS was sized to the peak heating lo- 
cation and the same thickness was used for the entire heat shield, resulting in a significant mass 
penalty. For the second, the TPS thickness was set piecewise given the expected heat profile 
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along the vehicle surface. Specifically, it was assumed that the heating fell off from maximum 
heating location to 60% of maximum value. Figure 9 shows how the TPS was divided into sec- 
tions with different thicknesses. Sections 2 and 3 are of constant thickness and sections 1 and 4 
are of variable thickness. Section 2 is 1.0 m long for the robotic case and 2.0 m long for the 
manned case. 

The third TPS sizing technique was a more aggressive version of thickness tailoring. In this 
case, no TPS was used on the back side of the sphere or cylinder sections. The maximum TPS 
thickness was the intersection of the sphere and cylinder at the centerline. The thickness decreas- 
es from maximum to 60% of the heating value 2.35 m down the cylinder. Along the circumferen- 
tial direction, the thickness decreases from maximum 60% of the heating value to the minimum 
thickness that can be manufactured, which is assumed to be 5.0 mm. Figure 10 shows how the 
heat shield is divided in this case. 



Figure 10. Aggressive TPS tailoring. 


RESULTS 

The overall trajectory design was an iterative process. First, the initial trajectory design space 
was generated and candidates were selected from the “pool” of initial conditions that satisfied the 
trajectory and mission requirements, e.g., presented low heat rates and g-loads. These candidates 
were deemed representative of nominal trajectories. These nominal trajectories, one for both the 
robotic and manned cases, were then used to size the TPS as previously discussed. The TPS siz- 
ing tool produces a new vehicle mass that is then used to re-run the trajectory design space, and 
the process is repeated until the TPS mass converges. Although further iterations would be re- 
quired to fully close the design, the current results shown here represent HAVOC trajectories that 
deliver the desired payloads to the target altitudes. 

Thermal Protection System 

Table 4 lists the materials and values for the constant thickness approach for the robotic mis- 
sion with a PICA heat shield at the end of the first iteration. Table 5 lists the materials and repre- 
sentative values for a HEEET heat shield for the same case. RTV-560 is a bonding agent. 


Table 4. Representative properties of a PICA heat shield, robotic mission 


Material 

Thickness (cm) 

Density (kg/m^) 

Mass (kg) 

CBE+25% (kg) 

PICA 

4.9 

273.9 

1651.4 

2640.3 

RTV-560 

0.03 

1409.6 

51.4 

64.2 

SIP 

0.229 

131.3 

36.5 

45.7 

RTV-560 

0.03 

1409.6 

51.4 

64.2 

Gr/BMI 

0.1 

1605.0 

194.9 

243.7 

A1 honeycomb 

5.08 

70.5 

434.8 

543.5 

Gr/BMI 

0.1 

1605.0 

194.9 

243.7 
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Table 5. Representative properties of a HEEET heat shield, robotic mission 


Material 

Thickness (cm) 

Density (kg/m^) 

Mass (kg) 

CBE+25% (kg) 

HEEET RL 

0.2 

1099.9 

271.4 

339.3 

HEEET IL 

2.79 

824.4 

2811.5 

3514.4 

RTV-560 

0.03 

1409.6 

51.4 

64.2 

Gr/BMI 

0.1 

1605.0 

194.9 

243.7 

A1 honeycomb 

5.08 

70.5 

434.8 

543.5 

Gr/BMI 

0.1 

1605.0 

194.9 

243.7 


For the PICA TPS mass (with margin) is 2,238.3 kg and the strueture is 1,030.9 kg, for a total 
of 3,269.2 kg. For the HEEET TPS mass (with margin) is 3,917.9 kg and the strueture is 1,030.9 
kg for a total of 4,948.8 kg. HEEET demonstrates a higher mass beeause it is more effieient at 
higher peak heat rates. Thus, for missions with relatively benign heating environments (sueh as 
the robotie HAVOC mission), in general PICA would be the preferred option due to its lower re- 
quired mass. 

A summary of the TPS sizing results using the three different teehniques is listed in Table 6. 
Both PICA and HEEET are viable eandidates for the robotie mission. PICA is the reeommended 
option due to its lower mass. In the ease of the manned mission, though PICA is lighter, it is ap- 
proaehing its heat flux limit. HEEET is also strueturally stronger than PICA, and ean tolerate 
mueh higher peak heat fluxes (on the order of 1,000 W/em^). Steepening the trajeetory would in- 
erease the peak heat flux but reduee the total heat load and reduee the overall mass of the HEEET 
system in both the robotie and manned mission seenarios. Thus, despite having a higher mass, 
HEEET is the reeommended option. The dual pulse eapability must be experimentally verified for 
both materials. 


Table 6. Summary of TPS masses. 



Robotic 

Manned 

Material 

Constant 

Moderate 

Aggressive 

Constant 

Moderate 

Aggressive 

PICA 

2,238.3 

2,124.6 

1,329.4 

14,151.2 

13,536.1 

11,813.2 

HEEET 

3,917.9 

3,527.3 

1,775.2 

30,645.8 

28,708.2 

25,066.0 


Using moderate TPS tailoring, the eurrent EV mass (payload, airship, and aeroshell), is ap- 
proximately 4,537 kg for the robotie mission (PICA) and 133.9 mt for the manned mission 
(HEEET). Using an aggressively tailored TPS permits a signifieant reduetion in the overall TPS 
mass, providing the option to either have additional available payload mass or reduee the overall 
size of the vehiele. 

Aerocapture Trajectory 

Aeroeapture inputs for the eurrent mass iteration are listed in the Table 7. The lift-to-drag ratio 
was ehosen to take full advantage of both the lift and drag area of the ellipsled vehiele, and the 
initial veloeity and target were ehosen based on the arrival opportunity from Earth and the desired 
initial EDI eonditions. 
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Table 7. Nominal inputs for each system. 



Robotic 

Human-Crewed 

Human-Cargo 

Lift-to-Drag Ratio 

0.52 

0.52 

0.52 

Initial Velocity 

11.3 km/s 

11.3 km/s 

11.3 km/s 

Target Orbit 

300 km circular 

300 km circular 

300 km circular 

Aeroshell Diameter 

4.7 m 

10m 

10m 

System Mass 

4,537 kg 

98 mt 

128 mt 


The nominal trajectories associated with each set of inputs produce g-loads and heat rates, 
shown in Figure 11, within the requirements of the system and the capability of the TPS material 
selected for this study. As the mass of the system increases from the robotic mission to the hu- 
man-crewed mission, the peak heat rate value increases to reflect the change in ballistic coeffi- 
cient. A heavier and larger diameter vehicle, entering the atmosphere at the same speed, requires 
more drag force in order to hit the same apoapsis target. To accumulate the necessary drag force, 
that vehicle must dive deeper into the higher density sections of the atmosphere, increasing the 
peak dynamic pressure and g-loads, and producing high peak heat rate values. Thus the peak heat 
rate values for the low mass robotic mission are much lower than those of the very high mass 
human missions. 




Relative Velocity (m/s) 


Figure 11. G-load and total heat rate values for each of the three systems. 

Each of these trajectories delivers the vehicle into a 300 km circular orbit, enabling the EDI 
portion of the mission. 

EDI Design Space 

Table 8 lists relevant parameters from valid trajectories using moderately tailored thermal 
protection systems. The bolded lines are those deemed representative of a HAVOC trajectory and 
were selected as nominal. Entry conditions are the independent parameters that were varied in the 
design space. Bank modulation values, designed to mitigate lofting and excessive aerodynamic 
loads, were also varied. Peak loads and parachute deploy conditions are obtained from the simu- 
lated trajectories (recall that parachute deploy conditions are defined to be at Mach 2.1). 
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Table 8. Sample valid unmanned and manned trajectories. 


Entry Conditions 

Bank Modu- 
lation 

Loads 

Parachute Deploy Conditions 
(Mach 2.1) 

Bank 

. of path 

^ attack angle 

Bank Bank 
(g) (alt) 

Peak Peak Total 

g- heat heat 

load rate load 

1 Dynamic Relative 

Altitude ^ 

pressure velocity 

deg deg deg 

deg deg 

g W/cm^ J/cm^ 

km Pa m/s 


Robotic 


0 

40 

-4.25 

0 

70 

6.105 

119.24 

9093 

79.27 

1494 

466.04 

0 

40 

-4.75 

0 

70 

7.490 

134.04 

8582 

79.29 

1486 

465.93 

0 

40 

-4.75 

20 

70 

7.731 

135.27 

8464 

79.34 

1470 

465.69 

0 

50 

-4.75 

20 

70 

8.332 

115.24 

6639 

82.46 

713 

452.35 

10 

50 

-4.25 

20 

70 

6.821 

104.86 

7093 

82.03 

791 

453.91 

10 

50 

-4.75 

20 

70 

8.366 

115.46 

6630 

82.47 

711 

452.30 

20 

50 

-4.75 

20 

70 

8.468 

116.13 

6606 

82.51 

704 

452.16 

20 

50 

-4.75 

20 

85 

8.468 

116.13 

5782 

79.17 

1529 

466.53 

20 

50 

-3.50 

20 

70 

4.845 

87.51 

7840 

80.13 

1230 

461.99 

Manned 

10 

50 

-4.25 

20 

70 

6.668 

273.65 

17760 

74.07 

4599 

486.80 


Several trends in the parameters of interest may be identified from the table (only one 
manned mission is shown, though the same trends seen in the robotie mission results hold for this 
elass of mission). Adjusting the entry angle of attaek affeets the final altitude: as angle of attaek 
deereases, the final altitude deereases due to the EV having less lift (reeall that peak lift is at a 55° 
total angle of attaek). Thus, higher angles of attaek are desirable to terminate at a higher altitude, 
where there is lower density and therefore lower dynamie pressure. 

Adjusting the entry flight path angle has a small effeet on final altitude, but a strong effeet on 
g-load: as the entry flight path angle approaehes zero, the g-load deereases, as expeeted. Thus, 
shallower entry flight path angles are desirable from a g-load mitigation standpoint, but steeper 
angles are desirable to avoid atmospherie skip-out. The entry flight path angle of -4.75 (robotie) 
and -4.25 (manned) listed in the table was seleeted beeause it was deemed a safe margin from 
skip-out while also providing survivable g-loads. 

Other trends may be identified from the table. Inereasing the entry bank angle inereases the 
peak g-load but has a negligible effeet on peak heat rate and final altitude. Inereasing the target 
bank angle for loft mitigation has a negligible effeet on peak heat g-load and peak heat rate, but 
deereases the final altitude, resulting in higher density and dynamie pressure. Inereasing the target 
bank angle for g-load mitigation has negligible effeet on peak heat rate and final altitude, but in- 
ereases the peak g-load. 

Thus, it is desirable to have an entry bank angle elose to vertieal, a high entry angle of attaek, 
a shallow flight path angle, a bank angle modulation for g-load mitigation elose to vertieal, and a 
bank angle modulation for loft mitigation elose to horizontal. This is supported by the design 
spaee results shown in Figure 6 and Figure 7. 

Finally, it ean be seen that none of the manned missions have a dynamie pressure at the Maeh 
2.1 eondition that is between 200 and 900 Pa, whieh is required for supersonie paraehute deploy. 
As will be diseussed later, this issue must be addressed with the use of an additional aerodynamie 
deeelerator. Nevertheless, it was of interest to examine these trajeetories in further detail. Analy- 
sis of the terminal deseent and airship inflation phase of the manned mission trajeetory eontinued 
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under the assumption that an additional deeelerator is used to establish feasible paraehute deploy 
eonditions. 

EDI Trajectory 

Trajeetory results for the nominal HAVOC missions are shown in Figure 12 through Figure 
14. Reeall that these are not intended to be the proposed trajeetory profiles, but rather show tra- 
jeetories drawn from the many valid eases that fulfill the mission requirements. Figure 12 shows 
trajeetory parameters for both elasses of HAVOC missions. Lofting behavior is drastieally miti- 
gated from initial trajeetories (though not eompletely eliminated) and brief altitude inereases are 
restrieted to below 100 km. 



Figure 13 shows aerodynamie parameters for both elasses of HAVOC missions. The most no- 
table differenees between the missions in these results are those between dynamie pressure and 
ballistie eoeffieient. The eause for the higher dynamie pressure for the manned mission is two- 
fold. Figure 12 showed that a given time, the manned mission is both at a lower altitude and mov- 
ing at a higher veloeity due to the higher ballistie eoeffieient. Reealling that = 2 dynam- 

ie pressure is sensitive to ehanges in veloeity. Combined with the higher density at lower alti- 
tudes, the dynamie pressure profile for the manned mission is signifieantly higher than that of the 
robotie mission. In faet, from Table 8, by the time the entry vehiele reaehes the nominal ehute 
deploy eondition of Maeh 2.1, the dynamie pressure is approximately 4.6 kPa. As will be dis- 
eussed later, this dynamie pressure is too high for paraehutes eommonly used in EDL operations. 

The two smaller peaks in dynamie pressure at approximately 6 and 7 minutes into the manned 
trajeetory are eaused by the slight lofting behavior seen in the top panel of Figure 12. 
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time (min) 


Figure 13. Aerodynamic parameters. 

Figure 14 shows aerodynamic and aeroheating loads for both classes of HAVOC missions. 
The higher heat rates and correspondingly higher heat loads for the manned mission are evident. 
G-loads, shown in the bottom panel, illustrate the effects of the bank angle modulations designed 
to mitigate these loads. 




time (min) 

Figure 14. Aerodynamic and aeroheating loads. 
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Recall that the dynamic pressure at Mach 2.1 for the case of the manned missions well ex- 
ceeded the 900 Pa limit of conventional supersonic parachutes. Despite this limitation, analysis of 
the terminal descent and airship inflation phase of the manned mission trajectory continued under 
the assumption that an additional decelerator may be used to establish feasible parachute deploy 
conditions. 

Airship-related parameters during terminal descent and airship inflation is shown in Figure 15. 
These are shown as a function of altitude rather than time. Note the logarithmic scales on volume 
and drag force. 





Drag force (N) 

Figure 15. Terminal descent and airship inflation parameters. 

The plots show several features of interest. Airship inflation begins when the descent speed 
becomes less than 100 m/s and then continues at a constant rate (listed in Table 2) until the air- 
ship is fully inflated. Full inflation occurs at an altitude different than the targeted 50 km. The 
neutral buoyancy condition of the airship, however, is designed to be 50 km, and so in both cases 
the airship will continue to descend or ascend until it settles at the target altitude. 

The drag force profiles exhibit two sudden decreases at 77.4 km and 67.5 km for the robotic 
case and at 62.6 km and 52.6 km for the manned case. These correspond to aeroshell and para- 
chute jettison, respectively. In the case of the robotic mission, the parachute is jettisoned just prior 
to full inflation (recall that the parachute is jettisoned when the buoyancy force reaches 90% of 
the parachute drag). Corresponding increases to the Mach number at these altitudes are also ap- 
parent. 

SUMMARY AND CONCLUSIONS 

A trajectory design and analysis that describes aerocapture, entry, descent, and inflation of 
manned and unmanned Venus lighter-than-air missions has been presented. Analysis of point 
cases drawn from the manned mission design space showed that dynamic pressures at parachute 
deploy exceed the valid environment due to the relatively low altitude at which the vehicle reach- 
es Mach 2.1, and resultant high density. The high dynamic pressure is beyond the survival capa- 
bility of commonly used parachutes (such as disc-gap-band) used in current planetary explora- 
tion. Thus, a different technology such a ribbon parachute or ballute must be used to slow the ve- 
hicle enough to permit the airship to inflate while under a parachute. In the case of the unmanned. 
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robotic lighter-than-air mission to Venus, it was determined that such a mission is feasible using 
current and near-future technology, such as a PICA thermal protection system. Future work will 
include implementation of guidance and control systems aboard the entry vehicle to aid in refin- 
ing the reference trajectories. For the manned mission this may improve the dynamic pressure at 
the Mach 2.1 condition and enable the use of a conventional supersonic parachute. 
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